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ABSTRACT 

We present deep optical 18-medium-band photometry from the Subaru telescope over the ~ 30' x 30' 
Extended Chandra Deep Field-South (ECDF-S), as part of the Multiwavelength Survey by Yale- 
Chile (MUSYC). This field has a wealth of ground- and space-based ancillary data, and contains 
the GOODS-South field and the Hubble Ultra Deep Field. We combine the Subaru imaging with 
existing UBVRIzJHK and Spitzer IRAC images to create a uniform catalog. Detecting sources 
in the MUSYC "BUi?" image we find ~ 40,000 galaxies with Rab < 25.3, the median 5a limit of 
the 18 medium bands. Photometric redshifts are determined using the EAZY code and compared to 
~ 2000 spectroscopic redshifts in this field. The medium band filters provide very accurate redshifts 
for the (bright) subset of galaxies with spectroscopic redshifts, particularly at 0.1 < z < 1.2 and at 
z > 3.5. For 0.1 < z < 1.2, we find a la scatter in Az/(1 + z) of 0.007, similar to results obtained 
with a similar filter set in the COSMOS field. As a demonstration of the data quality, we show that 
the red sequence and blue cloud can be cleanly identified in rest-frame color-magnitude diagrams at 
0.1 < z < 1.2. We find that ^20% of the red-sequence-galaxies show evidence of dust-emission at 
longer rest-frame wavelengths. The reduced images, photometric catalog, and photometric redshifts 
are provided through the public MUSYC website. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — astronom- 
ical databases: surveys — astronomical databases: catalogs 



1. INTRODUCTION 

Studies of distant galaxies and active galactic nuclei 
(AGN) require redshift measurements to determine dis- 
tances and associated lookback times. For optically- 
bright objects, redshifts can be measured effi ciently with 
multi-slit spectrographs on large telescopes ()Davis et al.l 
l2003tlLe Fevre et al.ll2005T) . However, spectroscopic red- 
shift measurements are very difficult for galaxies and 
AGN that are at z > 1, obscured, or intrinsically faint. 
As a result, we had to rely almost exclusively on pho- 
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tometric redshifts determined from broad- band photom- 
etry for galaxies with I > 23 (see, e.g., iGrazian et al.l 
20061: IWuvts et~all[2008l , and many other st udies) . 

The COMBO-17 survey (IWolf et alJ l200l pioneered 
the use of medium-bandwidth filters as a compromise 
between imaging and spectroscopy. These filters sample 
the spectral energy distributions (SEDs) of galaxies at a 
resolution of R — 10 — 20 and provide a redshift quality 
of 1-2%, intermediate between spectroscopy and broad- 
band imaging. This improvement enables measurements 
of rest-frame colors and the environment of galaxies, and 
improves the accuracy of determinations of the physical 
properties of galaxies. It also opens up the possibility to 
directly detect strong emission lines, which is particularly 
relevant for the identification of AGN. 

With the development of medium- band filter sets in 
the near-I R (Ivan Dokkum e t al. 2009|) and on 8m class 
telescopes (|Ilbert et al.ll2009f ). accurate redshifts are now 
becoming available for objects that are well beyond the 
limits of spectroscopy, and the full potential of this tech- 
nique is being realized. 

In this paper, we present deep 18-band optical 
medium-band photometry from the Subaru telescope in 
the - 30' x 30' "Extended" Chandra Deep Field-South, 
as part o f the Multiwavelengt h Survey by Yale-Chile 
fMUSYC: [Gawiser et al.ll2006bf l. Th e ECDF-S field ha s 
an extensive set of ancillary data (see iTavlor et al.12 009) , 
recently augmented by very deep Spitzer IRAC data in 
the SIMPLE^ survey (Damen et al. 2P1P, submitted). 
Furthermore, the field contains the GOODS-South field 
and the Hubble Ultra Deep Field. Our medium band sur- 

3 |http : //www. astroT yale . edu/dokkum/SIMPLE/ 
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vey is similar to that of Ilbert et al. (2009) in the larger 
COSMOS field, although we include 18 medium bands 
while there are only 12 available for the larger COSMOS 
field. 

We discuss the observations in Section [2] and the de- 
tails of data reduction in Section U3 We include ancil- 
lary data in our catalog from the MUSYC survey, which 
is described in Section [4] The details of our photome- 
try are discussed and the public catalog is presented in 
Section [5] We describe the photometric redshift deter- 
minations and evaluate the improvements provided by 
medium-band filters in Section [6l Finally in Section [7] 
we present the color-magnitude diagrams of the ECDF- 
S out to redshift 1.2 and investigate the occurrence of 
dusty galaxies on the red sequence. Throughout this pa- 
per we assume Hq = 71 kms -1 Mpc^ 1 , Ct m = 0.3 and 
n A = 0.7. 

2. OBSERVATIONS 

We carried out a medium-band imaging campaign on 
the Subaru Telesc ope using the wide field-of-view cam - 
era Suprime-Cam (|Mivazaki et al.ll2002t Hve et al.ll2004D . 
Observations were collected in 6 runs over 2 years, in Jan- 
uary, February and December 2006, and January, March 
and December 2007; the observations are summarized in 
Table [T] We show the response curves of the medium- 
band filters used in this study in Figure [TJ These filters 
were designed to maximize the accuracy of photometric 
redshifts by evenly sampling the wavelength range from 
400-900 nm with an effective resolution of A/AA ~ 23 
(|Taniguchill200l . 

Suprime-Cam consists of ten 2k x 4k CCD chips that 
cover a large field of view (24' x 27'). A single raw ex- 
posure in filter IA651 is shown in Figure [21 each of the 
10 individual chips show large gradients in background 
(from sensitivity variations and vignetting), as well as 
regions of bad pixels, and there are small gaps between 
the individual CCDs. Therefore, our observing strat- 
egy consisted of taking multiple Suprime-Cam exposures, 
dithering the pointing slightly to cover the small chip 
gaps. The median number of individual exposures for 
a given filter was 9, but the number of exposures and 
the exposure time were varied depending on conditions 
at the telescope and the available time each night. The 
long exposure on I A 738 was designed to be used as a 
detection image for relevant science projects, but for our 
purposes here we prefer the significantly deeper MUSYC 
BVR image as a detection image. Table Q] lists each fil- 
ter (1), the full- width at half maximum in nanometers 
(2), its observation date (3), the number of dithered ex- 
posures (4), and the total exposure time in seconds (5). 

3. DATA REDUCTION 

The data were reduced using a combination of stan- 
dard IRAF0 tasks; sdf red, a data reduction a nd analysis 
software package writte n for Suprime Cam ( Ya gi et al.1 
120021 lOuchi et al~ll2004h : and c ustom tasks. Our r educ- 
tion procedures foll ow those of lLabbe et al.l (|2003| ) and 
IQuadri et al.l (|2007|) . and are briefly described here. 

14 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under a cooperative agreement with 
the National Science Foundation. 



TABLE 1 
Medium Band Observations 



Band 


FWHM [nm] 


Date Observed 


# Exp Exptimc [s] 


(1) 


(2) 


(3) 


(4) ' (5) 



IA427 
IA445 
IA464 
IA484 
IA505 
IA527 
IA550 
IA574 
IA598 
IA624 
IA651 
IA679 
IA709 
IA738 
IA767 
IA797 
IA827 
IA856 
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Fig. 1. — Subaru medium-band filter transmission curves, includ- 
ing atmospheric transmission, quantum efficiency, and the trans- 
mission of the optical elements of the telescope plus instrument. 
This filter set was designed to sample the spectrum evenly be- 
tween 400 and 900 nm with an effective resolution of A/AA ~ 23 
in order to compute accurate photometric redshifts. 

3.1. Flat Fielding and Bias Subtraction 

Dome flats are well-exposed with high Signal-to-Noise 
ratios (S/N) and show the pixel-to-pixel variations on 
the CCD, while dark sky flats more accurately reproduce 
the variation in sensitivity to the night sky spectrum 
and the illumination pattern across the CCDs. There- 
fore, for each filter both dome flats and dark sky flats 
are constructed. We divide the high S/N dome flats by 
the dark sky flats, effectively removing the differences 
between the CCD response to the dome lights and the 
light of the night sky, i.e., flattening the dome flats. To 
maintain the high S/N, we then smooth the resulting 
flats with the iraf routine boxcar. The smoothing ker- 
nel for the boxcar routine was selected to be 10 pixels on 
a side because this minimizes the variance between ad- 
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Fig. 2. — Example raw data frame from Suprime-Cam, for a 
600s exposure in Filter IA651 taken on 2007-1-18 (East-up, North- 
right). Each exposure is composed of 10 CCDs in a 5 x 2 array 
with small gaps between each chip. There are large gradients in the 
background illumination of each chip from sensitivity variations 
and vignetting. There are bad pixel columns and regions of low 
sensitivity near the chip edges, which must be accounted for in the 
data reduction. 



jacent smoothed areas (<0.1%); this effectively increases 
our S/N by an order of magnitude. 

We also constructed master bias frames, combining 
^10-20 bias frames collected each observing run with the 
IRAF routine zerocombine. One master bias frame per 
run was sufficient as the bias frames were stable from 
night to night. 

After we construct the flat and bias frames, we apply 
them to the raw data frames. First, the two-dimensional 
master bias frames are subtracted from each raw data 
frame, then the remaining overscan correction is sub- 
tracted using sdf red. Next, flat fielding is done by first 
dividing the raw science exposures by the dome flat and 
then dividing this by our smoothed sky flat, as follows: 



raw frame — bias dome flat. 

v / ~Fi 4~ ) smoothed 



domeflat 



sky flat 



(1) 



Even after careful flat fielding, small gradients in the 
sky background can remain that can affect photome- 
try non-uniformly acros s the field. Using SExtractor 
(|Bertin fc Arnoutsl figQe). we measured and removed the 
sky gradient across the field using a large 2D mesh (1000 
pixels on a side). The large mesh was selected to avoid 
overestimating the sky near the edges of extended ob- 
jects and near the location of faint objects not detected 
in individual exposures. 

3.2. Image Combination 

After flat fielding and bias subtraction, we are left with 
multiple individual exposures for each filter, each com- 
posed of 10 chips, ready to be combined. These chips 
must be aligned and then combined, preserving the flux 
in the stars and rejecting bad pixel artifacts. 

Suprime-Cam provides a very large field-of-view, but 
with large geometric distortions caused by the optics. 
These distortions must be corrected for when combin- 
ing individual dithered data frames into a final image. 



For a first-order coordinate correction, we use the stan- 
dard task in sdfrd. Then we apply a secondary cor- 
rection using the deep BVR combined image from the 
MUSYC Survey (jGawiser et all l2006b| ) using the iraf 
routine mscimatcfcpl. 

In our image combination, we maximize the S/N in 
the seeing disk of point sources by combining the indi- 
vidual CCD chips from each exposure using a weighted 
average. The weighted average is determined by the 
RMS of the sky, the average FWHM and the relative 
flux scale for each of the individual CCD exposures, fol- 
lo wing the procedure des cribed in detail in Appendix A 
of iGawiser et all (|2006b| ). which was used in creating 
the broad-band point-source-optimized images for the 
MUSYC survey. This results in an optimal weight for 
each individual CCD chip as they are combined to make 
the final image. The weights are normalized to the first 
exposure central bottom chip. 

Additionally, to improve our image quality, we use 
carefully constructed bad pixel maps when combining the 
final images. Bad pixel maps are first created from pix- 
els marked by sdfred, including bad columns, isolated 
dead pixels and the CCD chip edges. Then careful visual 
inspection is made of each exposure and all bad pixels 
found by eye, in addition to artifacts such as satellite 
trails, are included in the final bad pixel maps. 

The final images are then created by averaging each 
of the individual weighted chips, rejecting the bad pixels 
and using the iraf routine combine with a percentile 
clipping algorithm. Experimentation showed the per- 
centile clipping algorithm provided optimal cosmic ray 
rejection. In Figure El we show final reduced images for 
I A 651. On the left of Figure we show the entire 
frame to show the flatness we achieved across the im- 
age, while the right panel zooms in on region displayed 
in the green box in the left panel. In the zoom-in, we see 
that many galaxies are clearly visible to the eye, the local 
background is relatively flat and no cosmic rays remain. 
We list details about each final combined image in Table 
d The FWHM, Column 2, is calculated using the iraf 
routine imexam on several hundred stars. 

3.3. Photometric Calibration 

We now have a single final image for each filter, 
but these must be placed on a standard flux scale. 
To achieve this, each night at the telescope a handful 
of ESO spectro-photometric standard star£3 were ob- 
served. These stars are used to calculate the zero-point 
for each final image. 

We calculate the magnitudes for these standard 
stars in our filter system {m ss ) by convolving the ESO 
standard star spectra with the effective transmission in 
each filter. We can then calculate the zero-point for 
each filter, using the amount of flux observed in each 
star (flux b s ) in the standard star images. The zero 
point for the final image in each filter is given by: 



— 2.5 x log flux bs — airmass x k v (2) 



15 The MSCRED package in IRAF was originally 
developed for mosaic reductions by NOAO. Informa- 
tion about the NOAO Mosaic Project can be found at 
http : //www.noao . edu/kpno/mosaic/mosaic . html 

10 http: / / www.eso.org/sci / observing/tools / standards /spectra/ 
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Fig. 3. — Left: Full frame showing the final combined image in filter IA 651, rotated 90 degrees from the raw frame shown in Figure[2] 
(North is up, East is to the left). Note the uniformity across the field from careful attention to flat fielding and the dithering pattern of 
individual exposures. The S/N at the edge of the field is lower due to the smaller number of individual exposures combined at the edges 
with our dither pattern. A box indicates the region of the field displayed in the right panel, zooming in by a factor of 5. Right: Detail of a 
final combined image in the filter IA651. Notice the detailed morphology visible for nearby large galaxies. 



Here fc„, the air-mass coefficient, stands for the extinc- 
tion in magnitudes per airmass and is necessary because 
each standard star observation was made at a differ- 
ent airmass from the science images. To determine fc„, 
we first obtained the Extinction Curve for Mauna Kea 
from the CFHT observers' manuaQ and compared this 
to an extinction curve der ived on Mauna Kea in the B 
and V bands published by IKrisciunas et aTl ([1987). Ad- 
ditionally, on nights when a single standard star was 
observed at multiple airmasses, we directly calculated 
extinction curves. All fc„ measurements were consis- 
tent with each other. Therefore when calculating the 
zero points for each filter, we used our own measured 
values for the six filters for which they were available 
(IA427, JA484, JA527, IAQ2A, IA738, an d IA8 56), two 
additional values from IKrisciunas et al.l (| 1987T) ; the B 
band for I A 445 and the V band for I A 550 and in- 
terpolated values for the remaining 10 filters according 
to the curve published by the CFHT. The fluxes are nor- 
malized such that they are fluxes per second, to account 
for the different exposure times of the standard stars and 
the final image exposures. The air-mass assigned to the 



final combined image is that of the first exposure, which 
is used to normalize the flux in the image combination 
( q3.2j) . The zero points for each filter are listed in Table 
El Column 4. 

3.4. Noise Properties 

The noise properties of each final image must be un- 
derstood in order to calculate the nominal depth of each 
final combined image, as well as accurate photometric 
errors. SExtractor assumes Poisson sky noise and thus 
can underestimate the total errors in the photometry, 
which include contributions from electron read-out noise, 
sky noise, imperfect background subtraction and pixel- 
to-pixel correlations that can be introduced during the 
reduction process. To accurately account for all of these 
effects, we sum the counts in a large number of aper- 
tures randomly placed throughout the image, and add 
these error estimates i n quadrature with the output er- 
ror from SExtractor (lLabbe et aTl 120031 : iGawiser et al.l 
l2006bt iQuadri et al.ll2007| ). The apertures are placed in 
locations to avoid objects (using SExtractor's segmenta- 
tion map) and have an identical size to the photometry 



apertures we later use in each band (See Section [5]). The 
http://ww.cfht.hawaii.edu/Instru^^ f fl ux j n these apertures 
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TABLE 2 
Medium Band Image Properties 



TABLE 3 



Band 


r WxlM [ J 


ocr acptn [ArSJ 


Zero Point [AB] 


(!) 


i 2 ) 


(3) 




(4) 


IA427 


1.01 


25.01 


25 


10 


± 0.11 


IA445 


1.23 


25.18 


25 


07 


± 0.08 


IA464 


1.79 


24.38 


25 


30 


± 0.03 


IA484 


0.76 


26.22 


25 


50 


± 0.05 


IA505 


0.94 


25.29 


25 


34 


± 0.02 


IA527 


0.83 


26.18 


25 


72 


± 0.03 


IA550 


1.13 


25.45 


25 


88 


± 0.06 


1I\0 i *± 


u.yo 


9^ 1 « 

iO. ID 


25 


71 


± 0.02 


IA598 


0.63 


26.05 


26 


02 


± 0.03 


IA624 


0.61 


25.91 


25 


89 


± 0.05 


IA651 


0.60 


26.14 


26 


15 


± 0.03 


IA679 


0.80 


26.02 


26 


20 


± 0.03 


IA709 


1.60 


24.52 


26 


02 


± 0.03 


IA738 


0.77 


25.93 


26 


02 


± 0.02 


IA767 


0.70 


24.92 


26 


04 


± 0.02 


IA797 


0.68 


24.69 


26 


02 


± 0.02 


IA827 


1.69 


23.60 


25 


92 


± 0.04 


IA856 


0.67 


24.41 


25 


73 


± 0.01 



a Total Magnitude 

describes the noise in the background on the image, and 
the sigma width of this Gaussian measures the uncer- 
tainty in the background noise (i.e., the depth) of the 
imaging. These depths are based on flux measurements 
within the aperture size used for photometry. Integrating 
over the PSF, an offset of —0.6 magnitudes is necessary 
to correct t hese magnitudes from aperture values to total 
values (see 35.2)1 . We report 5a depths with this correc- 
tion in TableEl Column 3. We note an independent anal- 
ysis of these same data give slightly different zeropoints 
and depths, due to small differences in methodology. The 
zeropoints are consistent to ~ 0.1 mag for all filters. 

3.5. Final Images 

We created final combined images for each filter (Fig- 
ure [3]), which are available from the MUSYC webpageP^I. 
We list each image and its properties in Table [2j Column 

1 lists the central wavelength of the filter in nm, Column 

2 lists the seeing, as measured using the average FWHM 
of a sample of several hundred stars in the actual image, 
Column 3 lists the 5 a depth in AB magnitudes, and 
Column 4 lists the zero point value for the image. 

4. ANCILLARY DATA 

4.1. MUSYC Survey: Existing Optical and 
Near- Infrared Data 



The Extended Chandra 
S) has been ta rgeted by 



frared surveys (lArnouts et al 



Deep Field-South (ECDF- 
host of optical and 



12001 IMov et al.l l200l 



Wolf et al.ll2004HGawiser et alj2006b[|Hildebrandt et al.l 
2006t Damen et al 2010, submitted). We o btain reduced 
optica l and near-infrared imaging from iTavlor et al.l 
(2009), where the observations, reductions and character- 
istics of these data are described in detail. To summarize 
briefly, the UU^BV RI imaging originate from the ESO 
archive and were combined from multiple projects us- 
ing the Wide Field Imager (WFI) on the ESO MPG 2.2 

18 http : //physics . rutgers . edu/~gawiser /MUSYC 
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a Total Magnitude 
1 Gawiscr et al. ( 2006. 



Ulildebr andt et aTl i PHOT 
^ Taylor eTall lt2009l 1 
IMov et alj H2003I) 
f Damen et al. 2010, submitted 

m telescope. These data were collected and calibrated 
as part of GaBoDS (th e Garching-Bonn Deep Survey; 
iHildebrandt et all I2006D . The z-band data, collected 
as part of the MUSYC survey, are from t he Mosaic-II 
came ra on the CTIO 4m Blanco telescope dMuller et al.l 
11998ft and are des c ribed further bv lGawiser et al.1 ()2006b) 
and ITavlor et al.1 J2009). The H- band data, taken wit h 
Soft on the ESO NTT 3.6 m telescope (IMov et al.| [200l, 
covers 80% of the field (jTavlor et al.1 12009TI . The JK 
imaging was obtained usi ng the ISPI came ra on the 
CTIO Blanco 4m telescope (jTavlor et alj|2009f ). The im- 
ages along with the size of the FWHM, 5tr depth, zero 
point and survey in which they were observed, are listed 
in Table [3J The FWHM and 5er depth are measured in 
the same way we determined these values for the Sub- 
aru data f ^3.4ll3"?2")) and are listed in Columns 2 and 3 
of Table [S) We adopt the ca libration and zero points 
published in ITavlor et al.ll2"009l shown here in Column 4. 

4.2. SIMPLE Survey 

Considerable Spitzer time has been invested in the 
ECDF-S with b oth IRAC (Infrared Array Camera; 
iFazio et al.l 120041 ) and MIPS (the Multi- Band Imaging 
Photometer for Spitzer, iRieke et al.ll2004T ). The Spitzer 
IRAC/MUSYC Public Legacy in the ECDF-S (SIMPLE, 
Damen et al. 2010, submitted ) project obtained very 
deep IRAC imaging across the full ECDF-S. The IRAC 
data cover 3 — 8 fim, to ~24th magnitude AB at 3.6^m. 
The astrometry was calibrated using the MUSYC BVR 
detection image, the same imag e we used in the Subaru 
medium-band reductions ( §3.2)1 : the resulting positional 
accuracy for individual sources is < 0.3 arcseconds (la; 
Damen et al. 2010, submitted). We include the IRAC 
data from the SIMPLE Survey in Table [3] The values 
for FWHM (Column 2) and 5 a depth (Column 4) are 
measured in the same manner as for the ground-based 
imaging ( ^3.41 [372]) . The zero point values (Column 4) 
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come from (Damen et al. 2010, submitted). 

4.3. Astrometric Calibration 

In order to perform multi-band photometry, each im- 
age is transformed to the image plane and pixel scale 
(0."2 67) of the stacked B VR image used for detec- 
tion (jGawiser et all 12006b). The RMS astrometric er- 
rors, compared to an astrometric catalog of the ECDF-S 
(Terry Girard, private communication), a re estimated to 
be less than 0.2" across the entire field (|Gawiser et al.l 
2006b). Using IRAF tasks geomap and geotran all images 
were resampled to the coordinate system of the BVR im- 
age (both in x-y pixel coordinates and in right ascension 
and declination) with a North- up tangent plane projec- 
tion. 

5. PHOTOMETRY 

In order to create accurate multiwavelength SEDs, a 
constant fraction of light needs to be collected for each 
object across every band. If the varying point spread 
functions (PSFs) of each band are not taken into account, 
similar apertures will collect different fractions of an ob- 
ject's total light. The data in our survey come from a va- 
riety of telescopes, with large variations in seeing (Tables 
[5] and El Column 2). Experiments with the broad band 
images showed that failure to correct for this effect would 
bias the U-V colors of sources by up to ~ 20%. Therefore, 
close attention was paid to obtaining accurate colors. In 
this section we describe the methodology used to mea- 
sure accurate colors f H5.ip . the selection of the apertures 
to be used for the photometry f ^5.2p . the detection of 
objects on the images ( H5.3I) and the completeness of our 
detection as a function of object magnitude ( ^5.41) . 

5.1. PSF Matching 

We have a total of 10 ground-based broad-band im- 
ages (U, C/38, B, V, R, I, z, J, H, K), 4 IRAC images 
(3.6 /im, 4.5 /im, 5.8 /im, 8.0 /im), and 18 medium-band 
images (IA427, IAU5, 7A464, JA484, 7A505, IA527, 
IA550, IA574, IA598, JA624, JA651, IA679, IA709, 
JA738, IA767, IA797, 7A856). The seeing in these im- 
ages ranges from ~ 0.5" to over 2" (Tables [5] and [3]). 
In this section, we describe our technique for smoothing 
multiple images to the same PSF to correctly measure 
color fluxes. 

One technique commonly used to achieve uniform pho- 
tometry across images in multiple filters relies on smooth- 
ing all images to the PSF size of the image with the 
worst seeing (e.g lLabbe et al.ll200l iCapak et all 120071; 
iTavlor et al.l 120091 ). Then large apertures are used for 
photometry, collecting flux within a uniform physical 
area across all images for each galaxy. However, degrad- 
ing all images to the largest PSF size decreases the S/N 
significantly in each filter where the image is smoothed. 
Given the large range in PSF size across our data, if we 
degraded all images to the worst seeing we would sacri- 
fice the excellent image quality found in the majority of 
our observations. 

Therefore, we apply a two-fold approach. First, for the 
12 images with narrow PSFs (Tables [H |3J, we smooth 
them to the PSF of the BVR ~ 0"8 image. For these 
images, all fluxes are measured using a single aperture, 
which in this case provides accurate PSF-matched pho- 
tometry. Then, for the images with PSFs larger than 



that of the BVR image (Tables HJ |3J , we degrade a copy 
of the BVR image to match each of the larger PSF sizes. 
We can then measure a color, e.g., f(K)-f(BVR,smoothedh 
and scale this to the aperture flux measured in the BVR 
0"79 image, fmvR)- I n other words, the PSF-matched 
K-band flux is measured by: 

, , f(BVR) , s 

J (K,psf matched) — J(K) x ~i \^) 

J (BV R. smoothed) 

Here fmvR) 1S the flux measured in the BVR image 
with its native 0"8 seeing, and fmvR,smoothed) is the nux 
measured in the copy of the BVR image, smoothed to 
have the same seeing as the if-band image. Therefore, 
f(K,psf matched) provides accurate colors when compared 
with the fluxes measured from the 12 images with nar- 
row PSFs. We note that the aperture radius in which the 
flux is measured is selected to be the FWHM (see tj5.2p . 
so to measure the flux in the iC-band image (fix)) an d 
the smoothed BVR image (fmv R,smoothed)) , we use an 
aperture radius equal to the FWHM of the iT-band im- 
age. We note that in Equation 3, K represents any of the 
filters whose final images have PSF sizes larger than the 
BVR image (U, U38, B, V, I, z, J, H, K, IAA27JAU5, 
IAA6A, IA505, IA550, IA57A, IA709, and 3.6, 4.5, 5.8, 
8.0 /im). The final catalog contains these PSF-matched 
fluxes for each of the 32 bands. 

Here we describe our technique for smoothing two im- 
ages to an identical PSF. We first attempted to smooth 
the images with a Gaussian convolution kernel, but found 
that the residual PSF variations were large due to the 
non-Gaussian PSF shape for stars in the images. In- 
stead, we decided to build a separate PSF for each filter 
directly from the images themselves. To do this, we se- 
lected ~ 20 isolated stars in each band that were well 
exposed but not saturated. The selected stars were reg- 
istered, normalized to their peak flux and then averaged 
to create a single model PSF for each band. We used the 
Lucy-Richardson algorithm (iRAF's lucy) to construct 
a kernel to convolve with each image. In Figure [4j we 
show the stellar curve of growth for the 12 images with 
PSFs smaller than the BVR image (dotted lines). The 
corresponding new curves of growth for the 12 smoothed 
images (red lines), where we have convolved the origi- 
nal image with our Lucy-Richardson kernel to match the 
BVR image PSF, show a consistent stellar profile. The 
ratio of the flux enclosed within our photometry aper- 
ture to total flux (measured within 6 FWHM) is stable 
to nearly ~ 1%. This accuracy is also achieved in the 
smoothed copy of the BVR image for each filter with 
poor seeing. 

5.2. Aperture Selection 

In order to optimize the S/N for the photometry, the 
flux in each object is measured in the central high surface 
brightness regions of the objects and then later corrected 
to a total flux measurement. For a Gaussian PSF and un- 
corrected noise, the aperture diameter t hat maximizes 
the S/ N is 1.35 times the seeing FWHM (jGawiser et al.l 
2006b]). We use an aperture radius equal to the seeing 
FWHM, which encloses ~50% of the flux (Figure gj) and 
provides nearly optim al S/N for point source photometry 
(jGawiser et a!]|2006b[ ). 

We note that an aperture defined to ensure high S/N 
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Fig. 4. — Stellar curve of growth for the 12 bands with PSFs 
smaller than the BVR PSF (black dashed lines), measured using 
20 isolated stars. Red solid lines show the smoothed PSF after con- 
volution with the appropriate kernel to match the BVR PSF, also 
indicated is the aperture used for the PSF-matched flux. Insert: 
BVR stellar curve of growth and fifth-order polynomial fit. We use 
this curve to measure the fraction of light missed in an aperture of 
a given size, and thus the correction to total flux ( i|5,3l l. 

can result in fluxes that measure only the central re- 
gion of galaxies, particularly for large galaxies at low 
redshifts. For example, at z=0.3, a 0'.'8 aperture is ~ 4 
Kpc. Therefore, for a Milky Way size galaxy at redshift 
0.3, our aperture would measure only the central bulge 
stars. Hence, any color gradients in the outer regions 
of the galaxy would not influence these aperture fluxes. 
Therefore it is important to note that in our catalog, the 
fluxes describe only the central regions of the galaxies 
at lower redshift. We also emphasize the importance of 
measuring the same central fraction of each galaxy in the 
apertures ( §5.1|) . 

When measuring the noise in the image (§ 13. 4|) . we 
used an aperture size equal to the seeing FWHM and 
converted these to total magnitudes assuming a stellar 
profile (Figure[4j. Therefore our 5c depths are computed 
for point sources, a magnitude offset of -0.6 from the 
aperture depth derived from the images. 

5.3. Source Detection 

For object detection w e run SExtractor version 2.4.4 
(jBertin &: Arnoutslll996[ ) in dual-image mode using the 
original BVR image for detection. In the 20 filters where 
a copy of the BVR image was smoothed to the image's 
larger PSF, SExtractor was run twice: first on the fil- 
ter's image and then again on the smoothed BVR image 
( M5. 1[) . SExtractor directly gives us aperture fluxes (with 
a radius of 2 x FWHM) for the images of the 12 filters 
with small PSFs. For the 20 filters with images with 
larger PSFs, we compute aperture fluxes from Equation 

m 

To correct our aperture fluxes to total fluxes for each 
source, we compute a total-flux correction factor for each 
object. SEx tractor's AUTO flux uses a flexible Kron-like 
(|Kronlll980l ) elliptical aperture to compute the flux. The 
flux in this aperture accounts for the size of the source 
and the source s hape, and measures appr oximately 94% 
of the total flux (jBertin fc Arnoutsl [l996;). We use SEx- 
tractor to obtain AUTO fluxes for each object from the 



BVR image. Then we correct this AUTO flux to a to- 
tal flux using the stellar curve of growth from the BVR 
image (Figure [4] insert, points). We fit this growth curve 
with a 5th order polynomial (Figure |4] insert, solid blue 
line) and compute the fraction of light enclosed for each 
object at the Kron radius (light frac). The total flux 
correction factor to be applied to the AUTO flux is then 
the inverse of the fraction of light enclosed at the AUTO 
flux radius for each source (totcor = 1/ light frac). To 
correct our aperture fluxes in any given filter to total 
flux in that filter we scale the aperture flux using the ra- 
tio of the BVR AUTO flux to the BVR aperture flux for 
that source and then multiply it by the flux correction 
factor (totcor): 

, , f BVR, AUTO , , 

JK,tot = J K, aperture X ■ X totcor , (4) 

J BV R. aperture 

where K can be replaced here by any filter in our catalog. 
We note that for extended obj ects this is a mini mum 
correction to a total magnitude (jTavlor et al.ll2009h . 

5.4. Completeness 

Although we measure the nominal depth of our images 
using the background noise ( fc|3.4[) in empty apertures 
across the field, the actual completeness of our catalog 
is a strong function of source magnitude. For each filter 
in the catalog, the 5cr depth is listed in Tables [2] and [3J 
for the BVR detection image, this depth is 26.8 in AB 
magnitudes. 

To empirically measure our overall effectiveness in de- 
tecting objects as a function of magnitude, we construct 
a stellar PSF, normalize it to various magnitudes, place it 
into the BVR detection image and determine if SExtrac- 
tor recovers the inserted star. The stellar image we use 
is that constructed from 20 well exposed stars in Section 
l5.Il to measure the BVR image PSF. A thousand random 
locations are selected, at locations at least 10 pixels (^3 
FWHM) away from detected objects and avoiding re- 
gions of the field washed out by bright stars. Near larger 
objects (eg. galaxies) , distances to these random sources 
were increased to be greater than twice the Kron radius. 
We create multiple images, each with 1000 stars of vary- 
ing magnitude and then run SExtractor on each of these 
new images. Figure [5] shows the percentage of these stars 
recovered as a function of simulated magnitude (solid 
line). The 95% completeness level is nearly 27 mag, and 
it falls quickly below that magnitude. This agrees well 
with the number determined from placing empty aper- 
tures around the field. However, in both cases we ignore 
the extent of detected objects on the field. Because real 
galaxies do overlap each other along the line of sight, we 
repeat this exercise, no longer restricting the random po- 
sitions to avoid detected objects of this field (although 
we still avoid the very bright stars which wash out large 
areas on the field). The recovery rate of these sources 
falls much faster at faint magnitudes, due to missing the 
inserted stars near brighter stars and galaxies (Figure [SJ 
dashed line). Still, at a magnitude of ~25.5 we are 90% 
complete in our detections. 

Figure [5] shows our BVR band number counts for 
the BVR-dected catalog. The number counts increase 
steadily until BVR~25.5 mag [AB], where we are still 
90% complete. Differential number counts can measure 
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Fig. 5. — The fraction of simulated stars detected by Sextractor. 
For stars placed to avoid bright stars and the locations of other 
objects on the image (solid line), the 90%/95% completeness level 
is nearly 27 AB magnitudes. Stars placed randomly (dashed line), 
yield lower completeness because they suffer from source confusion. 
The 90% (95%) completeness for these simulated stars is ~25.5 (24) 
AB magnitudes. 
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Fig. 6. — Number counts in the BVR-band. The number of 
sources increases steadily to BVR~ 25.5 mag, our 90% complete- 
ness level. Our differential number count measurement 0.339 ± 
0.004 consistent with the value measured by (Gawiscr ct al. 2006b) 
0.34 ± 0.01. 



the geometry of space and the evolution of structure in 
the universe. Our fit to the differential number counts 
per magnitude per square degree is 0.339 ± 0.004, con- 
sistent with the 0.34 ± 0.01 measured bv IGawiser et al.l 
(|2006bl ). 

5.5. Galactic Extinction 

The ECDF-S is at high Galactic latitude and there- 
fore has a very low Galactic extinction. For the location 
of the ECDF-S, RA=3 /l 32' & DEC=-27°48', we cal- 
culate a value of E(B - V) = 0.0088 from the 100 /xm 



maps of iSchlegel et al.l (11998 ^^1. We calculate the ex- 
pected Galactic extinction in each band, as suming R=3.1 
and us ing the Galactic Extinction Curve of Cardclli ct al. 
( 1989) with updates in the optical region from lO'Donnelll 
(1994). These values are listed in Table El Column 1. 
The Galactic extinction corrections are not included in 
the catalog, but are applied before photometric redshifts 
are computed (<j6]). 

5.6. Catalog 

We present our catalog format in Table 01 available 
online in full formatB The photometry is measured 
in units of flux (/xJy) and is not corrected for Galac- 
tic extinction. For the BVR combined image, the AUTO 
flux as well as the APERTURE flux is included. The 
fluxes are aperture fluxes and can be corrected to to- 
tal fluxes using the total flux correction and the AUTO 
flux (Equation [4]). The geometrical parameters are out- 
put from Sextractor including the Kron Radius, A_IMAGE, 
B_I MAGE, THETA_I MAGE and CLASS_STAR. The table is pre- 
sented both as a text file and as a FITS file. In the 
FITS version of the cat alog, SExtractor detection flags 
(jBertin fc Arnoutslll996[ ) for each filter are included. 

6. REDSHIFTS 

The Subaru imaging allows us to obtain accurate 
(Az/(1 + z) ~ 0.01) photometric redshifts for the sources 
in our catalog. In this section we describe our method 
for computing photometric redshifts and evaluate their 
accuracy as determined by a subset of sources with spec- 
troscopic redshifts. 

6.1. EAzY 

In order to obtain highly accurate photometric red- 
shifts, we used EAzY, a program optimized to pro- 
vide high quality redshifts over < z < 4, where com- 
plete spectroscopic cal ibration samples are not available 
i|Brammer et al.l l2008f ). EAzY is a full- featured red- 
shift fitting code, allowing for the use of priors in com - 
puting photometric redshifts (e.g., BPZ; Bcnftcz 1999). 
It includes a user-fried ly interface based on HYPERZ 
(jBolzonella et al.l HOOQ) and a carefully selected tem- 
plate set, designed to optimize photometric redshifts for 
optical-NIR surveys (|Brammer "eTall l2008h . The tem- 
plate set and the magnitude priors are based on semi- 
analytical models that are complete to very faint mag- 
nitudes, rather than highly biased spectroscopic sam- 
ples, and so are particularly useful for samples of ob- 
jects such as du st obscured galaxies that are faint in 
the optical (e.g., iTreister et aLll2009al ). When running 
EAzY, we allow for a linear combination of all templates 
and include a broad-line AGN template only when fit - 
ting sources detected in X-rays (jCardamone et alJ 2008) . 
The AGN template is based o n the SDSS QSO tem- 
plate from iVanden Berk et aTl (|2001l ). which we've ex- 
tended towards the re d using the mean QSO SEP from 
iRichards et all pOOl ). Following lllbert et al.l pOOl) . 
EAzY now includes strong narrow emission lines in its 
galaxy templates by estimating an [O II] emission line 
flux from the UV luminosity of the scaled template and 

19 http:/ /www. astro. princeton.edu/ schlegel/dust/data/data.html 

20 http : / /physics . rutgers . edu/~gawiser/MUSYC 
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TABLE 4 

Summary of Photometric Catalog Contents 



Column No. Column Title Description 



1 


num 


Squential Object Identifier, beginning from 


2,3 


ra, dec 


right ascension & declination (J2000; decimal degrees) 


4 


CLASS-STAR 


SExtractor parameter measuring stellarity of object 


5 


Kron radius 


Sextractor parameter measuring source size in a flexible aperture (pixels) 


6,7 


AJMAGE,B.IMAGE 


SExtractor parameter measuring major and minor axis of image profile (pixels) 


8 


THETAJMAGE 


SExtractor parameter measuring position angle measured counterclockwise from North (degrees) 


9 


totcor 


aperture correction to convert AUTO flux to total (fiJj) 


10,11 


f_auto_BVR,e_auto_BVR 


SExtractor AUTO flux and error (/iJy) 


12-76 


LX,e_X 


aperture flux and error in each filter (/iJy), including that measured for the BVR image 


77 


flag_X 


BVR SExtractor detection flag * 



*note: in the FITS version of the table, SExtractor detection flags are included for all filters. 



adopting fixed line ratios for [O HI/O II], J H/3/0 II], 
and [Ha/O II] as defined by IKennicuttl {1998). Further, 
EAzY introduces a template error function to account 
for wavelength-dependent template mismatch. 

EAzY provides multiple estimates of the photometric 
rcdshift, including z_peak, which we adopt in this work 
as the photometric redshift estimate. Because z_peak is 
marginalized over the full probability distribution, it can 
differ from the best estimate from a straight x 2 minimiza- 
tion when there are two widely-separated peaks in the 
rcdshift probability function by selecting the peak with 
the largest integrated probability. Additionally, it in- 
cludes an estimate of the quality of each photometric red- 
shift (Q z ) which combines the \ 2 of the template fit, the 
width of th e 68% confide nce interval and the BPZ odds 
parameter (|Bem tcz 1999) in such a way t hat Q z increases 
as an y of those parameters deteriorates ()Brammer et al.l 
2008). We recommend adopting a cut in Q z , when using 
the p hotometric redshifts for science (seelBramme r et al.1 
120081 and tj6.5[) , and in this work adopt the conservative 
requirement of Q z < 1 for all photometric redshifts used. 

When using EAZY to compute our photometric red- 
shifts, we include the 18 medium bands and the optical 
and near-infrared ground based coverage in addition to 
the IRAC data. Example SEDs are shown in Figure [7J 
The FWHM of the broad band filters is indicated by a 
red line. We compute photometric redshifts for over 99% 
of sources in our BVR-detected catalog using the default 
EAZY output, and adopting the default template error 
function and R-band photometric prior. For the less than 
1% of BVR-detected sources that lie on the imaging area 
in fewer than 5 filters, we do not compute photometric 
redshifts. 

6.2. Spectroscopic Redshifts 

We collect all available spectroscopic redshifts for 
sources in our catalog from the literature in order to 
quantify our photometric redshift accuracy. The quality 
of these redshifts varies widely, from sources with multi- 
ple spectral line measurements to sources showing only 
hints of a single spectral line. Matching these spectro- 
scopic redshifts to our catalog, there are almost ~4000 
unique objects with redshifts, ~1000 of which have mul- 
tiple published spectroscopic redshifts. For sources with 
multiple redshift determinations we select those with a 
higher quality flag, favoring redshift measurements with 



published quality flags over those without. In Table [5J 
we list all references from which unique spectroscopic 
redshifts were obtained, including the number of spec- 
troscopic redshifts used from the data-set, the source 
paper, the quality flags we used in determining the ac- 
curacy of our photometric redshifts and the median R- 
band magnitude (§ 16. 5j) . When evaluating the accu- 
racy of the photometric redshifts, we restrict the qual- 
ity of the spectroscopic redshifts used in the compar- 
ison. If we include spectroscopic redshifts with lower 
quality flags, it lowers the determined accuracy by fac- 
tors of 2 or 3, independent of rcdshift or magnitude, and 
it increases the fraction of outliers. Each of the redshift 
catalogs adopted a different selection tec hnique for the 
source s for which they obtained spectra: iBalestra et al.l 
(2010) selected galaxies at 1.8 < z < 3.5 for the VI- 
MOS Low Resolution Blue grism and galaxies at z < 1 
in addition to Lyman Break Galaxies at z > 3.5 in the 
Mediu m Resolution (MR) orange grism, iVanzella et al.1 
(2008) selected galaxies using color criteria and photo- 
metric redshifts at 0.5 < z < 2 and 3.5 < z < 6.3 
for VLT/FORS2 observations. iLe Fevre etall (|200l se- 
lected galaxies using a magnitude selection of Iab < 24, 
ICimatti et al.l (|2002D selected galaxies using a magnitude 
selection of K s < 20, Lira et al., in prep, selected galaxies 
using a variety of criteria including sources with X-ray- 
counterparts, Lyman Alp ha Emitters an d Lym an Break 
galaxies at2.7<z<3.6, iSzokolv et~all (l2004h sele cted 
galaxies with X-ray-counterparts, IKriek et al.l 12008) se- 
lected K-bright galaxies at z ~ 2.3 using the BzK 
and DRG color select i on (|Daddi et al.H200"4t iFranx et al.l 
120031 ). iTreister et"aH (l2009bD selected galaxies with X- 
rav-counterparts. lStrolger et al.1 (|2004D se lected high red- 
shift galaxies in a Supernova Search, iCristiani et al.l 
( 2000) selected galaxies using t he VLT UV- Visual Echelle 
Spectrograph to detect Lya, Ivan der Wei et ail (l200l 
selected galaxies at z ~ 1 from the CO MBO- 17 cata- 
log wi th a compact and regular shape, and ICroom et al.l 
(|2001[ ) selected compact objects with J-K colors redder 
than the stellar sequence. We note that many of these 
spectroscopic datasets have a larger number of total pub- 
lished redshifts; we include here only those which we have 
adopted as unique redshifts in our spectroscopic redshift 
catalog. Aditionally, other references contain spectro- 
scopic redshifts in the ECDF-S, but they are not men- 
tioned here unless we adopt at least one spectroscopic 
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Fig. 7. — Example SEDs with EAZY spectral fits overlayed. In black are the data points with la errors. The blue line shows the 
linear combination of template SEDs that best fit the data. In the left hand column, we illustrate two low rcdshift examples where the 
Balmer break is fit by the medium-band filters. In the right hand column, we illustrate two high redshift examples showing a break due to 
Lyman-alpha absorption. We highlight fits to faint sources, R ~ 24 in 3 of the 4 sources shown. Emission lines do influence the photometry. 



rcdshift into our unique listing. 

6.3. Zero Point Adjustments 

To perform SED fitting, we require highly accurate col- 
ors across the optical and near-infrared spectrum. Pho- 
tometric zero points (Tables [2] and El Column 4) were 
determined individually in each band through standard 
star measurements. However, small offsets of a few per- 
cent in adjacent bands can introduce significant color- 
offsets. This is a particular concern when a catalog, like 
ours, includes photometry from a variety of instruments 
taken over a period of many years. If we assume that 
apriori we know the spectral shapes of the galaxies and 
that they can be well fit by our synthetic SED templates, 
than any systematic offsets of observed fluxes in a given 
filter are due to a photometric zeropoint error in that fil- 
ter. Because our sources cover a wide range of redshifts, 
each filter samples a different region in every rest-frame 



galaxy SED and therefore any systematic photometric 
offsets in this filter are due to zero point offsets rather 
than a template error at a single wavelength. 

To cross-calibrate our photometry, we take the best- 
fit EAzY SED for each galaxy and measure the offset of 
observed flux from that template flux in each filter. For 
normally distributed errors in flux measurements, the av- 
erage of these offsets should be zero. However, we find 
systematic flux offsets between the observed flux and the 
template flux for each filter. We wish to correct the fil- 
ters zero point by the median offset in each filter. For 
this calculation we restrict our sample to sources with 
spectroscopic redshifts and high S/N detections, in this 
case expected SED template fluxes of at least 5 x the 
detection limit for the filter. We calculate the median 
flux offset for each filter, iterating by applying the cal- 
culated flux offsets and re-fitting the SEDs. Within 3 
iterations, these flux offsets vary by less than 1%, and 
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TABLE 5 

Spectroscopy 



# Sources References 




(li ml ^'i~\r PI a ctq tyioH mn i R *-ni 


(1) (2) 




\ a ) w 


1239 Balcstra ct al. 2010 


A 23 


573 Vanzella et al. 2008 


A,B 24 


223 Le Fevre et al. 2004 


4,3 23 


224 Cimatti et al. 2002 


1,0 23 


211 Lira et al., 


in prep 




52 Szokolv et al. 2004 


3,2,1 23 


9 Kriek et al. 2008 


n/a 24 


7 Treister et al. 2009b 


1 23 


5 Stroleer et al. 2004 


n/a 25 


3 Cristiani et al. 2000 


n/a 24.5 


3 van der Wei et al. 2004 


n/a 26 


2 Croom et al. 2001 


n/a 22 




TABLE 6 




Photometric Zero Point Offsets from SED Fitting 


Band Galactic Extinction 


ZP Correction 


(1) 


(2) 


(4) 


U38 


0.041 


-0.184 


U 


0.043 


-0.245 


B 


0.034 


0.053 


V 


0.028 


0.075 


R 


0.022 


0.075 


1 


0.014 


0.032 


z 


0.013 


0.156 


.1 


0.008 


-0.023 


H 


0.005 


-0.011 


K 


0.003 


0.255 


IA427 


0.037 


0.285 


IA445 


0.036 


0.124 


IA464 


0.034 


-0.057 


IA484 


0.032 


0.075 


IA505 


0.030 


0.146 


IA527 


0.029 


0.010 


IA550 


0.027 


0.126 


IA574 


0.026 


0.054 


IA598 


0.024 


0.064 


IA624 


0.023 


0.116 


IA651 


0.022 


0.032 


IA679 


0.021 


-0.035 


IA709 


0.020 


-0.001 


IA738 


0.019 


0.054 


IA767 


0.018 


0.086 


IA797 


0.017 


0.075 


IA827 


0.016 


-0.106 


IA856 


0.014 


0.106 



The region around the MUSYC-Extended Chandra 
Deep Field- South has been observed by many surveys 
and therefore there are m any public catalogs available 
in this region. We refer to iTavlor et al. (2009) for a de- 
tailed list of these surveys and a quantitative comparison 
of their photometry. Here we briefly compare our cata- 
log to two prev i ous ver sions of the MUSYC cata log, from 
iGawiser et all (|2006bf ) and ITavlor et all (|2009D . and to 
the upda ted photometry o f the Combo- 17 medium-band 
catalog, (jWolf et al.ll2008f ). Each of these catalogs used 
slightly different techniques to derive the PSF-matched 
photometry, but all used Sextractor for detection and the 
same publicly available data . On average our fluxe s are 
within 1% of those quoted bylGawiser et all (l2006bl). For 



the co mparison to ITavlor et al.l ( 20091) and iWolf et al.l 



lf2008h . we include the offsets calculated in 36.31 Our 
photo metry agrees with that published by ITavlor et al.1 
(2009) on aver age to ~ . 02 ma gnitudes, and with that 
of Combo-17 (IWolf et al.l l2008h with an offset of -0.05 
magnitudes. 

6.5. Photometric Redshifts 

Comparing non- X-ray sources with high quality 
spectroscopic redshifts, we find a median accuracy 
Az/(1 + z) - O.OOcO out to z~5 (Figure [i left). Limit- 
ing ourselves to the subsample of sources with photomet- 
ric redshifts 0.1 < z < 1.2, where the Balmer break falls 
within our medium bands, we find Az/(1 + z) ~ 0.006 
(Figure [H right). We report the NMAD F^l f or cqmpar- 
ison with the works of others (e.g.. Illbert et ahl 120061 : 
iBrammer et~aLll2008h ITavlor et al.ll2009f ). and refer to it 
with the symbol tr z ~ 0.01. Wc define outliers as 10xcr z , 
and find that ~5% of the overall sample are labeled as 
outliers. 

In Table [Jj we show the photometric redshift quality 
as a function of galaxy magnitude; listing the magnitude 
bin, number of objects found within the given magnitude 
bin, the percent of the galaxies in the bin with Q z < 1 , 
the median 68% confidence interval for the galaxies in 
the bin as determined by EAzY, the number of galax- 
ies in the bin with high quality spectroscopic redshifts, 
the median Az/(1 + z), er z and the percent of galaxies in 
the bin with spectroscopically confirmed redshifts that 
are outliers. The number of high quality comparison 
spectroscopic redshifts is a strong function of magni- 
tude, but overall we maintain median Az/(1 + z) < 0.01 
(°z < 0.02) to an R-band magnitude of 25. The scatter 
in the photometric vs. spectroscopic redshifts almost cer- 
tainly irnder^s_tonates the ac tual uncertainty in the red- 
shifts Bra mmer et all (|2008| ). and therefore we use the 
median 68% confidence intervals computed by EAzY to 
compare the quality of the photometric redshifts for all 



we find that slightly altering the list of sources used to 
calculate these flux offsets does not affect our solutions 
by more than 1-2%. Therefore our photometry is sys- 
tematically uncertain at this level. For each band, we 
apply these photometric offsets, listed in Table El before 
computing final photometric redshifts used in this paper. 
We do not calculate offsets for the IRAC bands where the 
templates have greater uncertainties. 

6.4. Flux Comparison to other Catalogs 



Az — \z s pec Zpfootl 

22 Wc calculate the normalized median absolute deviation 



(NMAD) as 



"NMAD 



1.48 X median 



Az — median(Az) 



1 



as in (Brammer ct al. 2008). The normalization factor of 1.48 en- 
sures the NMAD of a Gaussian distribution is equal to its standard 
deviation and the subtraction of median(Az) corrects any offset 
from zero of the Gaussian. The NMAD is a useful measure of 
dispersion because it is less sensitive to outliers th at the standard 
deviation Ullbert et al.l|2006l ; iBTammer et aT1l2008T) 
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TABLE 7 

Photometric Redshift Quality vs Source Magnitude 



Magnitude 


# Objects 


Qz < l 


168 z — u68 z 


jf %spec 


Az/(1 + z) 


CT Z 


outliers 


17 < R < 19 


283 


99% 


0.011 


20 


0.004 


0.005 


5% 


19 < R < 20 


385 


98% 


0.012 


32 


0.005 


0.006 


3% 


20 < R < 21 


733 


98% 


0.014 


81 


0.006 


0.008 


1% 


21 < R < 22 


1676 


97% 


0.015 


167 


0.006 


0.007 


2% 


22 < R < 23 


3537 


93% 


0.017 


374 


0.006 


0.008 


2% 


23 < R < 24 


8087 


86% 


0.022 


524 


0.007 


0.011 


6% 


24 < R < 25 


17048 


68% 


0.038 


415 


0.014 


0.020 


5% 


25 < R < 26 


25628 


31% 


0.063 


77 


0.019 


0.032 


6% 


26 < R 


26127 


6% 


0.121 


7 


0.021 


0.026 


14% 


X-ray Sources 


825 


69% 


0.022 


236 


0.008 


0.012 


12% 



TABLE 8 

Photometric Redshift Accuracy vs Redshift 



Zphot 




# Qz < l 


168 z - u68 z 


Qz < 1 (% of Zspec) 


■fr Zspec 


Az/(1 +z) 


fz 


outliers 


o < 


z 


31381 


0.034 


92% 


1697 


0.008 


0.011 


4% 


o < 


z < 1.2 


22318 


0.024 


98% 


1242 


0.006 


0.008 


3% 


1.2 < 


z < 3.7 


8547 


0.087 


80% 


433 


0.019 


0.027 


8% 


3.7 < 


z 


516 


0.065 


92% 


21 


0.012 


0.016 


5% 



galaxies. This confidence interval is the difference be- 
tween the 68% confidence upper and lower bounds on 
the photometric r edshift computed from the probability 
distribution p(z) (|Brammer et al.1 [20081 ) and is a strong 
function of galaxy apparent magnitude. We find good 
quali ty photometric redshifts (Q z < 1; iBrammer et al.1 
120081 ) down to R-band magnitude of 25 for 70% of the 
sample (90% of the sample has Q z < 3), but miss a 
large fraction of the galaxy population at fainter magni- 
tudes. For many of the remaining sources poor fits are 
obtained by EAzY because the photometry is too un- 
certain, due to faintness of the sources, or intrinsic vari- 
ability in the source ov er the time period of which the 
photometry was taken (jSalvato et al.l I2009T ) . There are 
also cases where the intrinsic SED may not be matched 
by the templates and/or degeneracies in color-z space re- 
sult in multiple peaks in the redshift-probability distribu- 
tion. We note that when considering X-ray counterparts 
(Figure [3 purple diamonds), our accuracy is maintained 
Az/(1 + z) ~ 0.009, but the outlier fraction increases to 
18%. We include the X-ray counterparts in the last line 
of Table [7] for comparison. 

Additionally, we compare the photometric redshift ac- 
curacy as a function of redshift in Table [SJ Because the 
spectroscopic sample contains a variety of selection ef- 
fects ( N6.2p . Table [5] includes the median 68% confidence 
intervals computed by EAzY for all objects with Q z < 
1. We use photometric redshift to select the sources in 
each redshift bin so that we are comparing comparable 
sources in our measurements of 168 z — u68 z and cr z . Com- 
paring to the galaxies with spectroscopically determined 
redshifts, the completeness of the photometric redshift 
determinations (Q z < 1) varies as a function of redshift. 
We fit high quality photometric redshifts to over 90% 
of galaxies, but we note that the completeness falls to 
80% in the interval between 1.2 < z < 3.7. Overall, with 
medium-band photometry we have achieved highly accu- 



rate photometric redshifts for the majority of the sources 
in the ECDF-S. The photometric redshifts are available 
online through the MUSYC webpagj^l. 

We quantify the effect of the medium-band filters, com- 
paring the photometric redshifts with and without the 
additional information these filters provide. When ex- 
cluding the medium-band filters, we retain the broad 
band optical filters, as well as the J, H, K band and 
IRAC photometry. Overall the medium-band filters pro- 
vide a factor of three improvement over photometric red- 
shifts using broad-band filters alone, decreasing the me- 
dian Az/(1 + z) from 0.026 to 0.008 (comparing sources 
with Q z < 1). We note that with the addition of the 
medium band filters, high quality ( Q z < 1) photomet- 
ric redshifts are obtained for 20% more of the spectro- 
scopic redshift sample. The improvement is most no- 
ticeable at 0.1 < z < 1.2 where the median|Az/(l + z)| 
falls from 0.022 to 0.006 and at z > 3.7 where the me- 
dian Az/(l + z) falls from 0.024 to 0.006. These red- 
shift intervals are where the Balmer break (3700 A) and 
Lyman limit (912 A) fall inside the region of the spec- 
trum covered by the medium-band filters. However, we 
still see improvement at 1.2 < z < 3.7, where the me- 
dian Az/(1 + z) falls from 0.032 to 0.016. In Figure El 
we compare the spectroscopic sources with good fits to 
the SEDs (Q z < 1) with (left) and without (right) the 
medium-band filters. The medium-band filters not only 
tighten the accuracy around the z_phot=z_spec line, but 
also can help to rule out false redshift solutions (so called 
catastrophic failures) for sources with z_spec > 1.2. In 
each of these cases at 2 ~ 2.5, the broad band photom- 
etry alone is fit with by a Balmer break feature in the 
SED at low redshift (z ~ 0.1) but when the medium band 
photometry is added, the same optical region of the SED 

23 http : //physics . rutgers . edu/~gawiser/MUSYC 
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spec spec 

Fig. 8. — Left: Comparison of photometric and spectroscopic redshifts for all high quality spectroscopic redshifts (Tabled and photometric 
redshift fits with Q z < 1. The galaxies are colored by R-band magnitude, and those detected in X-rays are also indicated by purple diamonds. 
We find Az/(1 + z) ~ 0.008, or Az/(1 + z) ~ 0.006 at 0.1 <z<1.2, where the Balmer break falls into the wavelength range covered by our 
medium-band filters. Dotted lines are set at lOxNMAD, our definition of outliers. Right Zoom-in on low redshift region with linear scale. 



is fit by a Lyman alpha decrement caused by absorption 
by the IGM and the photometric redshift increases to a 
value more consistent with the spectroscopic determina- 
tion. An example is shown in Figure I1Q[ for a source 



with Rab ~ 24 and z s 



2.5. 



6.6. Star / Galaxy Separation 

We use two methods to identify the stars in our catalog , 
the first using a Bz'K color selection (1 Taylor et al.ll2009f) 
and t he second fitting stellar SED templates ( Ilbert et al.l 
l2009h . 

The Bz'K diagram is traditionally used to select mod - 
erate redshift (z > 1.4) galaxies (jDaddi et al.1 120041 ) . 
but is also an ef ficient discriminat i ng between stars 
and galaxies (e.g. iDaddi et all 12004 iBlanc et al.l 120081: 
iTavlor et al.ll2009l) . Beca use our filter s et is slightly dif- 
ferent than that used by IDaddi et al.1 (120041) . we apply 
the offsets determined bv IBlanc et al.l (|2008f ) of -0.04 
mag in (z-K) and 0.56 mag in (B-z) to our colors be- 
fore plotting them in the BzK diagram. In Figure [TTJ 
we show the BzK diagram for at all sources in our cat- 
alog with K AB < 21.84 (i.e., K vega < 20). The stellar 
sequence is clearly separated from the galaxies by the 
colors z — K < 0.3(5 — z) — 0.5, shown as a solid black 
line in Figure [TTJ For comparison, we include the hand- 
ful of spectroscopically identified stars (large blue + in 
Figure ITT1) . Only 2 of the 108 spectroscopically identified 
(Kab < 21.84) stars are missed by the BzK star color 
selection. Although this method of color selection is ef- 
fective at identifying stars, we are limited to the sources 
bright enough to be detected in our relatively shallower 
K-band coverage. 

To identify stars in all of our objects with photometric 
redshifts, we evaluate a x 2 value for z = stellar SED 
template (Pickles 1998) fits for each object. Since the 



IPicklesi ()1998[ ) templates do not go past 2.5/im, we do 
not include the Spitzer IRAC bands in the stellar fits. 
We flag as potential stars objects with reduced xltar — 
Xgaiaxy These sources are included in the BzK diagram 
as blue crosses. The stars identified through template 
fitting agrees well with those selected in our BzK color- 
selection, ~90% of stellar-template stars with Kab < 
21.84 have BzK colors of stars (Figure [TTJ dark blue 
crosses). Additionally, 80% of .Bzif-color-selected stars 
are also identified as stars by the stellar template fit. 
Finally, we look at the subset of sour ces detected in the 
GEMS imaging (Haussler et al. 2007). Selecting sources 
with stellar FWHM from the GEMS imaging, -80% of 
them are selected by our template fitting method as stars, 
similar to the 84% of point like sources that t he COSMOS 
team identifies as stars through SED fitting (jllbert et al.l 
2009). Furthermore, looking at bright (R < 24) extended 
sources in the GEMS imaging, we find fewer that 2% are 
misidcntified as stars by this template fit classification. 

6.7. Comparison to Combo- 17 Photometric Redshifts 

Combo- 17 produced the first highly accurate photo- 
metric redshifts of the ECDF-S region, and we wish 
to compare our own photometric redshifts to theirs 
(jWolf et al.ll2004D . We restrict ourselves to the brighter 
photometry (r < 24) where their published photomet- 
ric redshifts are more reliable (10% accuracy), restrict- 
ing ourseleves to r <22, we have fewer sources for com- 
parison, but the accuracy of Combo-17 (2%) approaches 
our own (Wolf et al. 20Q3). To examine our consistency, 
we compare our photometric redshifts to those deter- 
mined in Combo-17 for all sources contained in both 
catalogs, finding Az/(l + z) ~ 0.011 at R < 22 and 
Az/(l + z) ~ 0.029 at R <24. Therefore, our photo- 
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Fig. 9. — Comparison of photometric redshift verses z_spec for all sources with high quality spectroscopic redshifts and photometric 
redshift fits with Q z < 1. Note that the left panel, which includes the additional medium band filters in the fit, lowers the overall dispersion 
between the spectroscopic and photometric redshifts. 
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Fig. 10. — An example source at z~2.5 and Rab ~24, red lines indicate the FWHM of the broad band filters. Here the medium-band 
photometry increases the photometric redshift, ruling out the erroneous low-redshift solution found by fitting the SED with only broad-band 
photometry. In the case of the broad-band photometry alone, a Balmer break is fit to the UBV bands; when the medium band photometry 
is added, this same region is found to fit the sharper Lyman alpha decrement instead. 
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Fig. 11. — Bzk color-color diagram indicating stars. Note in 
this figure the z-K magnitude has been adjusted -0.04 and the B-2 
magnitude by + 0.56 to match the filter set used bv lDaddi et al.l 
(2004). Sources with spectroscopic redshifts 2 = (turquoise + ) 
and those which fit stellar SEDS (dark blue X), mostly lie along the 
stellar sequence in the bottom region on the figure. Also indicated 
are high redshift BzK galaxies including star forming sBzK galaxies 
(top left) and passive pBzK galaxies (top right). Overall, stars 
selected by template fitting to stellar SEDs are consistent with 
those indicated here on the BzK diagram. 

metric redshift determinations are consistent with those 
determined by Combo-17 to within their errors. 

To compare our respective accuracies, we use spectro- 
scopic sources (Figure [12]), and compare z_spec with pho- 
tometric redshifts for Combo-17 and our MUSYC survey. 
For R < 22, we find Az/(l + z) - 0.011 for Combo- 
17 and Az/(l + z) - 0.005 for MUSYC, and consider- 
ing all sources R A b <24 we find Az/(l + z) ~ 0.025 
for Combo-17 and Az/(1 + z) ~ 0.007 for our MUSYC 
redshifts. Therefore we find that although our photo- 
metric redshifts agree with those from Combo-17 sur- 
vey to within their respective uncertainties, the larger 
number of deeper medium-band filter observations in the 
MUSYC survey improves the photometric redshift accu- 
racy and extend this accuracy to fainter source popu- 
lations. We note that these differences come not only 
from the deeper medium-band data, but also by includ- 
ing JHK and IRAC photometry. 

7. THE RED SEQUENCE TOZ-l 

Here we demonstrate the quality of our data by show- 
ing the color-magnitude relation as a function of red- 
shift in the E-CDFS. The colors and morphologies of 
galaxies form a bim odal population dBlanton et alj[2003t 
IBaldrv et al.l [2004) . This bimodality represents a fun- 
damental relationship between mass and star formation 
history in galaxies, because the color of a galaxy is a 
proxy for the stellar age, where blue galaxies are younger 
and red galaxies are older. The red sequence contains 
the brightest galaxies, which typically have early- type 
morphologies while the blue cloud galaxies typically have 




z 

spec 



Fig. 12. — Comparison of our photometric redshifts (blue dia- 
monds) to those of Combo-17 (black X). Where our photometric 
redshift is closer to that determined from spectroscopy than that 
from Combo-17, a green line connects the two photometric redshift 
determinations. Where the photometric redshift from Combol7 is 
closer to that determined from spectroscopy, a red line connects the 
two photometric redshift determinations. Overall our photometric 
redshifts Az/(1 + z) are more accurate for many sources that were 
outliers in Combo-17. 



late- type morphologies. Relatively few galaxies are ob- 
served in the green valley, suggesting that galaxies spend 
only a short time here. How the bimodal galaxy popu- 
lation was created is one of the outstanding questions in 
galaxy studies today. 

Rest- frame fluxes were com puted following the proce- 
dure of iBrammer et afl (|2009T ). which measur es the rest- 
frame fluxes from the best-fitting template ()Wolf et al.l 
2003). This procedure is very different from a K- 
correction as it uses the observed medium bands that 
are closest in observed wavelength to the redshifted rest- 
frame band of interest. The Brammer et al. (2009) 
methodology is embedded in the EAZY photometric red- 
shift code (Brammer et al. 2008). 

In Figure [TBI we show the distribution of galaxies in 
the rest-frame U — V versus V plane. This plot most di- 
rectly relates observables to the bi-modal color sequence 
as a function of redshift. The restframe U — V color cov- 
ers the Balmer break and therefore is a good measure of 
recent star formation and the luminosity in the V band 
is a proxy for stellar mass. In F igure [T3l w e also include 
the red sequence cut defined bv lBell et al.l ()2004l ) in the 
Combo-17 survey (green dashed line) and our complete- 
ness limits in each redshift bin (blue solid line). We see 
a sequence of red points separated from the main blue 
clump of galaxies at all redshifts in our sample, high- 
lighting the quality of our data. The locati on of this 
red se quence is consistent with that found bv lBell et al.l 

(pool . 

7.1. Passive Galaxy Evolution 
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Fig. 13. — Rest-frame color-magnitude diagrams for galaxies as a function of redshift . The reds hift bins represent equal steps of co-moving 
volume, doubled in the last 3 bins. The green line is the red sequence cut defined by Bell ct al. ( 2004) and the dotted line shows the 90% 
completeness limit. The red sequence is clearly present out to z~1.2, where our redshift accuracy begins to fall. The red sequence contains 
the most luminous galaxies, while the blue cloud dominates the overall number density, especially for the fainter galaxies. 
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In color space both passive galaxies and dusty galax- 
ies can appear red in U — V colors. To investigate 
the nature of the red sequence (Figure [13j) , we use the 
IRAC photometry to separate out red passively evolv- 
ing galaxies from red, dusty, star-forming galaxies. The 
rest-frame V — J color can distinguish between these 
two populations because dust-free galaxies are blue in 
V- J color, while the dust-obscured galax ies are still 
red (lLabbe et all 120051: iWuvts et all |2007[) . Following 
IWilliams et"alT (|2009). we use a color-color diagram, rest- 
frame U-V vs V-J, to investigate the nature of the bi- 
modal galaxy color sequence for all galaxies with 0.2 < 
z < 1.2 (Figure . We color-code the galaxies from 
the r ed sequence on t he color-magnitude diagram red 
using [Bell et al.l (|2004l) 's red sequence cut. As expected, 
we find that passive and star f orming galaxies separ ate 
cleanly in color-color space (see, Williams et al. 20091) . 

We find that just over 20% of the galaxies identified 
as red sequence members in Figure Q2] have rest-frame 
V — J colors which place them in the star-forming se- 
quence, consistent with previous work at these redshifts 
using morphological informa tion (e.g.. IRuhland et aLl 
2009) and with the results in (jBrammer et al.ll2009D and 
Whitaker et al. (2010, in prep). Splitting our redshift 
interval in two (0.4 < z < 0.9 and 0.9 < z < 1.2) the 
fraction of red sequence galaxies with dusty colors is con- 
sistent between the two bins (21% and 23% respectively). 

8. SUMMARY 

In this paper we present new deep 18-medium-band 
photometry in the well-studied E-CDFS field. We 
reduced the raw data for the medium-band Subaru 
Suprime-Cam observations using a combination of stan- 
dard routines and custom tasks. The (public) catalog 
includes photometry from 10 ground-based broad-band 
images (17, U38, B, V, R, 7, z, J, H, K), 4 IRAC images 
(3.6 /zm, 4.5 /urn, 5.8 /im, 8.0 /im), and 18 medium-band 
images (IA427, 7^4445, IAAQA, 7A484, 7A505, 7,4527, 
7A550, 7A574, 7A598, 7A624, I Ami, 7A679, 7A709, 
7A738, 7A767, 7A797, IA856). The full catalog pro- 
vides multiwavelength SEDs for ~80,000 galaxies in the 
ECDF-S down to R [AB] ~ 27 (40,000 at R<25.2, the me- 



dian depth of the medium band imaging), although the 
accuracy and completeness of the photometric redshifts 
declines at R > 25.5. 

We computed accurate photometric redshifts using 
EaZy, a public photometric redshift code. The addi- 
tion of the medium band proves a factor of four improve- 
ment in the photometric redshift accuracy over the use 
of broad band filters alone. Comparing to spectroscopic 
redshifts, we find a scatter in Az/ (1 + z) of 0.008 for the 
full sample, 0.006 at 0.2 < z < 1.2 where the Balmcr 
break is covered by the medium band filters and 0.01 at 
z > 3.7 where the Lyman Limit (912A) is covered by the 
medium band filters. We find that the additional filters 
in the optical region even improve the photometric red- 
shifts for sources at 1.2 > z > 3.7, decreasing the median 
Az/(1 + z) from 0.03 to 0.02. 

We demonstrate that these photometric redshifts are 
sufficiently accurate to determine precise rest-frame col- 
ors of the galaxies. We detect the bimodal galaxy distri- 
bution out to z ~ 1.2 and find that 20% of the galaxies 
on the red-sequence have longer wavelength colors con- 
sistent with being dusty. In keeping with the spirit of 
the MUSYC collaboration, we provide public access to 
the images, photometry, and photometric redshifts. In 
future papers we will combine these accurate redshifts 
with the extensive ancillary data in the ECDF-S. 
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Fig. 14. — Rest-frame color-color diagram, showing U — V vs V — J (Labbe et al. 2005, Williams et al. 2009). A bimodal color-sequence 
is clearly visible separated by the dark lines ((Williams et al. 2009). Dust-free quiescent gala xies are blue in V — J, therefore the peak of 
galaxies at red U — V and bluer V — J are passive red sequence galaxies. All galaxies with the Bell et al. (2004) red sequence cut are shown 
as darker points. This diagram identifies 20% of those red sequence galaxies as dusty star forming galaxies rather than containing passive 
older stellar populations. 
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